Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation.

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

Disrupted chromatin architecture in olfactory sensory neurons: looking for the link from COVID-19 infection to anosmia.

We tackle here genomic mechanisms of a rapid onset and recovery from anosmia-a potential diagnostic indicator for early-stage COVID-19 infection. Based on previous observations on how olfactory receptor (OR) gene expression is regulated via chromatin structure in mice, we hypothesized that the disruption of the OR gene expression and, respectively, deficiency of the OR function can be caused by chromatin reorganization taking place upon SARS-CoV-2 infection. We obtained chromatin ensemble reconstructions from COVID-19 patients and control samples using our original computational framework for the whole-genome 3D chromatin ensemble reconstruction. Specifically, we used megabase-scale structural units and effective interactions between them obtained in the Markov State modelling of the Hi-C contact network as an unput in the stochastic embedding procedure of the whole-genome 3D chromatin ensemble reconstruction. We have also developed here a new procedure for analyzing fine structural hierarchy with (sub)TAD-size units in local chromatin regions, which we apply here to parts of chromosomes containing OR genes and corresponding regulatory elements. We observed structural modifications in COVID-19 patients on different levels of chromatin organization, from the alteration of whole genome structure and chromosomal intermingling to reorganization of contacts between chromatin loops at the level of topologically associating domains. While complementary data on known regulatory elements point to potential pathology-associated changes within the overall picture of chromatin alterations, further investigation using additional epigenetic factors mapped on 3D reconstructions with improved resolution will be required for better understanding of anosmia caused by SARS-CoV-2 infection.

AlloMAPS 2: allosteric fingerprints of the AlphaFold and Pfam-trRosetta predicted structures for engineering and design.

AlloMAPS 2 is an update of the Allosteric Mutation Analysis and Polymorphism of Signalling database, which contains data on allosteric communication obtained for predicted structures in the AlphaFold database (AFDB) and trRosetta-predicted Pfam domains. The data update contains Allosteric Signalling Maps (ASMs) and Allosteric Probing Maps (APMs) quantifying allosteric effects of mutations and of small probe binding, respectively. To ensure quality of the ASMs and APMs, we performed careful and accurate selection of protein sets containing high-quality predicted structures in both databases for each organism/structure, and the data is available for browsing and download. The data for remaining structures are available for download and should be used at user's discretion and responsibility. We believe these massive data can facilitate both diagnostics and drug design within the precision medicine paradigm. Specifically, it can be instrumental in the analysis of allosteric effects of pathological and rescue mutations, providing starting points for fragment-based design of allosteric effectors. The exhaustive character of allosteric signalling and probing fingerprints will be also useful in future developments of corresponding machine learning applications. The database is freely available at: http://allomaps.bii.a-star.edu.sg.

Conservation and Diversity in Allosteric Fingerprints of Proteins for Evolutionary-inspired Engineering and Design.

Hand-in-hand work of physics and evolution delivered protein universe with diversity of forms, sizes, and functions. Pervasiveness and advantageous traits of allostery made it an important component of the protein function regulation, calling for thorough investigation of its structural determinants and evolution. Learning directly from nature, we explored here allosteric communication in several major folds and repeat proteins, including α/ß and ß-barrels, ß-propellers, Ig-like fold, ankyrin and α/ß leucine-rich repeat proteins, which provide structural platforms for many different enzymatic and signalling functions. We obtained a picture of conserved allosteric communication characteristic in different fold types, modifications of the structure-driven signalling patterns via sequence-determined divergence to specific functions, as well as emergence and potential diversification of allosteric regulation in multi-domain proteins and oligomeric assemblies. Our observations will be instrumental in facilitating the engineering and de novo design of proteins with allosterically regulated functions, including development of therapeutic biologics. In particular, results described here may guide the identification of the optimal structural platforms (e.g. fold type, size, and oligomerization states) and the types of diversifications/perturbations, such as mutations, effector binding, and order-disorder transition. The tunable allosteric linkage across distant regions can be used as a pivotal component in the design/engineering of modular biological systems beyond the traditional scaffolding function.

Allosteric perspective on the mutability and druggability of the SARS-CoV-2 Spike protein.

Recent developments in the SARS-CoV-2 pandemic point to its inevitable transformation into an endemic disease, urging both refinement of diagnostics for emerging variants of concern (VOCs) and design of variant-specific drugs in addition to vaccine adjustments. Exploring the structure and dynamics of the SARS-CoV-2 Spike protein, we argue that the high-mutability characteristic of RNA viruses coupled with the remarkable flexibility and dynamics of viral proteins result in a substantial involvement of allosteric mechanisms. While allosteric effects of mutations should be considered in predictions and diagnostics of new VOCs, allosteric drugs advantageously avoid escape mutations via non-competitive inhibition originating from alternative distal locations. The exhaustive allosteric signaling and probing maps presented herein provide a comprehensive picture of allostery in the spike protein, making it possible to locate potential mutations that could work as new VOC "drivers" and to determine binding patches that may be targeted by newly developed allosteric drugs.

Exploring the Allosteric Territory of Protein Function.

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

Three-dimensional chromatin ensemble reconstruction via stochastic embedding.

We propose a comprehensive method for reconstructing the whole-genome chromatin ensemble from the Hi-C data. The procedure starts from Markov state modeling (MSM), delineating the structural hierarchy of chromatin organization with partitioning and effective interactions archetypal for corresponding levels of hierarchy. The stochastic embedding procedure introduced in this work provides the 3D ensemble reconstruction, using effective interactions obtained by the MSM as the input. As a result, we obtain the structural ensemble of a genome, allowing one to model the functional and the cell-type variability in the chromatin structure. The whole-genome reconstructions performed on the human B lymphoblastoid (GM12878) and lung fibroblast (IMR90) Hi-C data unravel distinctions in their morphologies and in the spatial arrangement of intermingling chromosomal territories, paving the way to studies of chromatin dynamics, developmental changes, and conformational transitions taking place in normal cells and during potential pathological developments.

AlloSigMA 2: paving the way to designing allosteric effectors and to exploring allosteric effects of mutations.

The AlloSigMA 2 server provides an interactive platform for exploring the allosteric signaling caused by ligand binding and/or mutations, for analyzing the allosteric effects of mutations and for detecting potential cancer drivers and pathogenic nsSNPs. It can also be used for searching latent allosteric sites and for computationally designing allosteric effectors for these sites with required agonist/antagonist activity. The server is based on the implementation of the Structure-Based Statistical Mechanical Model of Allostery (SBSMMA), which allows one to evaluate the allosteric free energy as a result of the perturbation at per-residue resolution. The Allosteric Signaling Map (ASM) providing a comprehensive residue-by-residue allosteric control over the protein activity can be obtained for any structure of interest. The Allosteric Probing Map (APM), in turn, allows one to perform the fragment-based-like computational design experiment aimed at finding leads for potential allosteric effectors. The server can be instrumental in elucidating of allosteric mechanisms and actions of allosteric mutations, and in the efforts on design of new elements of allosteric control. The server is freely available at: http://allosigma.bii.a-star.edu.sg.

AlloMAPS: allosteric mutation analysis and polymorphism of signaling database.

AlloMAPS database provides data on the causality and energetics of allosteric communication obtained with the structure-based statistical mechanical model of allostery (SBSMMA). The database contains data on allosteric signaling in three sets of proteins and protein chains: (i) 46 proteins with comprehensively annotated functional and allosteric sites; (ii) 1908 protein chains from PDBselect set of chains with low (<25%) sequence identity; (iii) 33 proteins with more than 50 known pathological SNPs in each molecule. In addition to energetics of allosteric signaling between known functional and regulatory sites, allosteric modulation caused by the binding to these sites, by SNPs, and by mutations designated by the user can be explored. Allosteric Signaling Maps (ASMs), which are produced via the exhaustive computational scanning for stabilizing and destabilizing mutations and for the modulation range caused by the sequence position are available for each protein/protein chain in the database. We propose to use this database for evaluating the effects of allosteric signaling in the search for latent regulatory sites and in the design of allosteric sites and effectors. The database is freely available at: http://allomaps.bii.a-star.edu.sg.

Exploring chromatin hierarchical organization via Markov State Modelling.

We propose a new computational method for exploring chromatin structural organization based on Markov State Modelling of Hi-C data represented as an interaction network between genomic loci. A Markov process describes the random walk of a traveling probe in the corresponding energy landscape, mimicking the motion of a biomolecule involved in chromatin function. By studying the metastability of the associated Markov State Model upon annealing, the hierarchical structure of individual chromosomes is observed, and corresponding set of structural partitions is identified at each level of hierarchy. Then, the notion of effective interaction between partitions is derived, delineating the overall topology and architecture of chromosomes. Mapping epigenetic data on the graphs of intra-chromosomal effective interactions helps in understanding how chromosome organization facilitates its function. A sketch of whole-genome interactions obtained from the analysis of 539 partitions from all 23 chromosomes, complemented by distributions of gene expression regulators and epigenetic factors, sheds light on the structure-function relationships in chromatin, delineating chromosomal territories, as well as structural partitions analogous to topologically associating domains and active / passive epigenomic compartments. In addition to the overall genome architecture shown by effective interactions, the affinity between partitions of different chromosomes was analyzed as an indicator of the degree of association between partitions in functionally relevant genomic interactions. The overall static picture of whole-genome interactions obtained with the method presented in this work provides a foundation for chromatin structural reconstruction, for the modelling of chromatin dynamics, and for exploring the regulation of genome function. The algorithms used in this study are implemented in a freely available Python package ChromaWalker (https://bitbucket.org/ZhenWahTan/chromawalker).

AlloSigMA: allosteric signaling and mutation analysis server.

MOTIVATION: Allostery is an omnipresent mechanism of the function modulation in proteins via either effector binding or mutations in the exosites. Despite the growing number of online servers and databases devoted to prediction/classification of allosteric sites and their characteristics, there is a lack of resources for an efficient and quick estimation of the causality and energetics of allosteric communication. RESULTS: The AlloSigMA server implements a unique approach on the basis of the recently introduced structure-based statistical mechanical models of allosteric signaling. It provides an interactive framework for estimating the allosteric free energy as a result of the ligand(s) binding, mutation(s) and their combinations. Latent regulatory exosites and allosteric effect of mutations can be detected and explored, facilitating the research efforts in protein engineering and allosteric drug design. AVAILABILITY AND IMPLEMENTATION: The AlloSigMA server is freely available at http://allosigma.bii.a-star.edu.sg/home/. CONTACT: igorb@bii.a-star.edu.sg.

First-principles study of heat transport properties of graphene nanoribbons.

We use density-functional theory and the nonequilibrium Green's function method as well as phonon dispersion calculations to study the thermal conductance of graphene nanoribbons with armchair and zigzag edges, with and without hydrogen passivation. We find that low-frequency phonon bands of the zigzag ribbons are more dispersive than those of the armchair ribbons and that this difference accounts for the anisotropy in the thermal conductance of graphene nanoribbons. Comparing our results with data on large-area graphene, edge effects are shown to contribute to thermal conductance, enhance the anisotropy in thermal conductance of graphene nanoribbons, and increase thermal conductance per unit width. The edges with and without hydrogen passivation modify the atomic structure and ultimately influence the phonon thermal transport differently for the two ribbon types.